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Abstract 

In this paper we consider the appUcation of electromagnetic theory to 
the analysis of the Scanning Near-field Optical Microscope (SNOM) in 
order to predict experimentally observable quantities such as the trans- 
mission or reflection coefficients for a particular tip-surface configuration. 
In particular we present the first application of a transfer matrix based 
calculation to this challenging problem by using an adaptive co-ordinate 
transformation to accurately model the shape of the SNOM tip. 

We also investigate the possibility of increasing the transmitted light 
through the SNOM tip by introducing a metal wire into the centre of 
the tip. This converts the tip into a co-axial cable. We show that, in 
principle, this can dramatically improve the transmission characteristics 
without having a detrimental effect on the resolution. 



1 Introduction 

One of the first lessons in optics which is drummed into our heads is that the 
resolution of any optical instrument is limited, fundamentally, by the wave- 
length of the light it uses. However, this resolution limit is not an absolute and, 
in fact, applies only to the image in the far field. Sub-wavelength information 
can, in principle, be obtained by probing the near field which contains both the 
propagating modes and evanescent modes bound to the surface of the object. It 
is these evanescent modes which correspond to sub- wavelength spatial frequen- 
cies and so encode the nano-metric detail of the object. So the principle which 
lies behind near-field imaging techniques is this: find some way to probe the 
near-field and imaging details smaller than the wavelength becomes feasible. 

In the early nineties, Pohl and co-workers at IBM [Q, |^, designed a near- 
field 'optical stethoscope' for imaging at a wavelength of 448nm with a claimed 
resolution of about A/20. Since then, many other groups have built Scanning 
Near-field Optical Microscopes (SNOM, or sometimes NSOM) along similar 
lines. A typical, modern SNOM can be described as follows. The surface under 
investigation is illuminated from above by laser light passing down an optical 
fibre. The end of the fibre is formed into a taper and the sides of the taper are 
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coated in about lOOnm of a metal, typically aluminium or silver the very end 
of the tip being left uncoated. In this way apertures can be reliably produced 
with diameters of ^ 20-500nm. As the tip is scanned over the surface, a shear- 
force mechanism is used to regulate the tip surface distance [Q, ||] . At each point 
the light transmitted through the sample is collected by a photomultiplier and a 
complete image built up. Such a device also offers the possibility for performing 
local, high resolution spectroscopy ^ || . 

Many variations on this basic idea are possible and there are several compre- 
hensive reviews . Possibilities include measuring the reflected rather than 
the transmitted light, illuminating from below and using the coated fibre tip as 
a light collector, or even illuminating from below above the critical angle and 
using an uncoated tip as a nano-scatterer to disrupt the near-field and frustrate 
total internal reflection |l2) . 

However, if the SNOM is to mature as an analytical tool then theorists have 
a critical role to play, both in order to interpret correctly the images produced 
by the microscope, but also to assist in the optimisation of tip and illumination 
system. Several approaches are possible for solving the problem of the elec- 
tromagnetic scattering between tip and surface, amongst those suggested have 
been finite difference time domain schemes and boundary matching methods. 
The former usually require a super-computer to handle a problem as complex 
as this. The latter, which expands the fields in each region in some suitable 
complete basis set and matches the fields at the boundaries, faces difficulties 
when handling regions of irregular shapes, which are needed for a realistic de- 
scription of the tip. However, one group have completely dominated the theory 
of the near-field microscope by introducing an iterative Green's function based 
method ||l3|, |l^. Using this method Girard etal. have been able to simulate 
images produced by reflection microscopes scanning over a range of dielectric 
and metal objects. 

2 A novel approach 

However, despite the sophistication of Girard's method it does rely on a real 
space discretisation in order to represent both object and tip. For an accurate 
description of the tip this is can be a serious restriction, especially if we are 
restricted to a regular, Cartesian grid. The problem is this: in reality the surface 
of the tip is a cone or paraboloid shape with a smooth surface, discretisation 
leads to a surface covered in steps, the size of which depends on the mesh size. 
The presence of these steps will result in additional scattering from the surface 
of the tip and may reduce the accuracy of the calculation especially if the size 
of the discretisation mesh is large. The only way to reduce the error is to use 
more mesh points leading to larger, slower calculations. 

We can, however, tackle the SNOM problem from a slightly different an- 
gle which helps us to overcome these problems. The idea is to perform the 
discretization in a co-ordinate system which follows the shape of the tip itself. 
While, in principle, the construction of such an adaptive, discrete mesh may 
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seem a horribly complex task, in practise this is not so. This is because of a 
useful result which allows us to factorise the problem into two parts. This re- 
sult, which is derived in detail elsewhere |]l5|, allows us to transfer all the details 
of the mesh geometry into an effective, tensorial permittivity and permeability 
(£«^, /t"^). Now we need only perform the discretization on a regular, cartesian 
mesh and then map that onto the adaptive mesh by introducing the appropriate 
£ and jl. 

How do Maxwell's equations appear in the adaptive co-ordinate system? The 
result is surprisingly simple. Let (71,(72,93 be the co-ordinates in our adaptive 
co-ordinate system and let the unit vectors Ui , U2 , U3 point along the generalised 
91j92,(73 axes. In this system Maxwell's equations become: 
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Notice that in general u^, and g"^ as well as e and /i will be functions of 
position. 



3 Details of our calculations 

We wish to calculate the transmitted intensity through a given tip and surface 
configuration and to this end we shall employ a transfer matrix method familiar 
to the theory of photonic band gap materials [^. Put simply this involves 
integrating the electric and magnetic fields from one side of the system to the 
other by means of a transfer matrix which relates the fields in one layer of the 
mesh to the fields in the next layer. 

As a first attempt calculation we will try to simplify the problem as much as 
possible. By exploiting the obvious cylindrical symmetry of the tip we can re- 
duce the problem from a three dimensional to a two dimensional one. The other 
major simplification we can impose is in the boundary conditions. By imposing 
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cylindrical symmetry we have confined the calculation to a cylindrical region of 
space containing both tip and surface. As usual for a transfer matrix method 
we can send an incident wave into the system from either end of the cylinder 
and find the transmission and refiection coefficients but we must specify 
the boundary conditions for the fields at the cyhnder's edge (see figure]^). The 
simplest form these boundary conditions can take is if we enclose the cylinder 
with a perfect metal, in other words we set the fields to zero at the boundary. 
This is not a trivial approximation to make, a much better solution would be to 
expand the fields outside the scattering region in some complete basis set and 
match the fields at the boundary, but it is not a foolish approximation either. 
As long as the wavefield does not have a sizable magnitude at the boundary and 
remains tightly localised between the tip and sample, as we might expect it to, 
then the choice of boundary conditions should not have a major effect on the 
results. 



3.1 An analytic form for the tip 



Our next choice is that of the discretisation mesh itself for this will determine the 
shape of the tip which we model. Ideally, we wish to specify the mesh geometry 
by means of an analytic function - this will greatly simplify the process of 
defining the effective permittivity and permeability, e and /t, which will contain 
all the information about the co-ordinate system. We also need to include 
sufficient parameters in the chosen function to allow us to fine tune the shape 
of the tip. 

One possible form for the mesh is this: 
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where ZR, 18 and IZ are the discrete variables which label each point on the 
mesh. This form contains five free parameters. The first three a,b,c describe 
an undistorted mesh. That is to say if (5 = a and a = 1 the mesh defined is 
just discrete cylindrical polar co-ordinates with a radial spacing of a, angular 
spacing of b and a spacing c in the z-direction. The other two parameters a 
and 6 control the distortion of this cylindrical mesh into a kind of hour glass 
shape; 6 defines the minimum radial distance between points at the centre of 
the mesh and a controls how tightly the other co-ordinate lines bunch in at the 
centre. Figure || shows the shape of the co-ordinate lines for a typical choice of 
parameters. 

The tip itself can now be defined by assigning the dielectric constant at each 
point in the mesh. The tip and substrate are glass, e ~ (2.16,0.0) the tip is 
coated with a 90nm layer of aluminium e = 1 — uj^/uj{ijJ + ij) with ujp — IS.leV 
and 7 = 0.27eV. The very end of the tip is not coated and this forms the 
aperture. A layer of metal around the top of the tip prevents stray light from 
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the source leaking around the side of the tip. Any object we wish to image 
can be introduced by changing the dielectric constant just above the surface. 
Figure ^ shows how the metal coating layer, substrate and object fit onto the 
distorted co-ordinate system. 

4 Results 

We will want to make comparisons between different tip designs to see if we can 
suggest ways to improve the utility of the near-field microscope. In order to do 
this we shall need some way to estimate the resolution of a given microscope 
configuration as well as to calculate the transmitted or reflected signal. One way 
to test the resolution is to compare the signals for two similar objects. Ideally, 
we would like to calculate the transmission for the tip in arbitrary positions 
above an arbitrary object and in this way build up a complete image of the 
object. Unfortunately, we have constrained ourselves somewhat by insisting 
upon cylindrical symmetry so this is not possible. The only choices available to 
us for the objects under the tip are rings and discs with the tip sitting directly 
above their centre. This is made clear in figure |[ Fortunately this is enough 
to get an estimate of the resolution, as we can compare a ring and a disc of the 
same outer radius to see whether we can detect the presence of the hole at the 
centre of the ring. 

Figures ^ to ^ show comparisons between transmission coefficients for ring 
and disc objects for three different tip apertures, given a fixed incident wave 
of frequency 2.5eV. The objects are made of the same glass as the substrate, 
54nm high and the tip-object separation is fixed at 36nm. Each graph gives the 
transmission against the outer radius of the object. The inner radius of the ring 
is a few nanometres smaller than the outer radius. Strictly, the outer radius is 
one radial mesh point larger than the inner radius so, because the mesh points 
are more closely packed at the centre, the difference increases as the size of the 
object increases. 

As the aperture decreases the overall magnitude of the transmission coef- 
ficient is reduced, from around 10~^ for an aperture radius of 48nm down to 
about 6 X 10~^ for a 25nm aperture. Reducing the aperture further proved diffi- 
cult as the magnitude of the transmission coefficient is reduced so much that the 
variations in it become comparable with the rounding errors in the calculation. 

4.1 Improving the design: Coaxial tips 

In order to get a high resolution, clearly what is required is as small an aperture 
as possible. However, because all the guided modes of the tip run into cut-off 
once the diameter is smaller than the wavelength, small apertures come at the 
price of a very small transmitted signal. This is borne out by the results in the 
previous section, where the transmission for an aperture of 25nm is about two 
orders of magnitude down on the transmission for a 48nm aperture. 



5 



However there is, in principle, a way to overcome this problem. Inserting a 
wire into the centre of the tip would convert the metal clad waveguide into a 
tapered coaxial cable. The coax is well known to have a different kind of guided 
mode called the TEM mode arising from the different boundary conditions 
satisfied at the centre. This mode has a free space like dispersion with no cut- 
off and so will propagate down the cable for any value of its radius offering the 
possibility of delivering a much higher photon flux to the end of the tip. 

There are one or two potential problems with using a coaxial tip in a mi- 
croscope, of course. The main one is the fear that the extra absorption caused 
by introducing a lossy metal into the centre of the guide would negate any gain 
given by the new mode, especially as the field distribution in the TEM mode 
behaves like 1/r so is largest along the central wire. We do our best to min- 
imise this by limiting the central wire only to the very end of the tip, where 
the radius has fallen below the cut-off point. We also try to maintain a uniform 
impedance along the coaxial part by ensuring that the central wire tapers in 
the same proportions as the rest of the tip. The other constraint we must be 
careful with is not to make the wire too thin. Once the radius of the wire is 
much less than the skin depth it will no longer be effective. 

In order to optimise the transmitted power through the coaxial tip we sought 
a compromise between the lack of cut-ofF in the TEM mode and the increased 
absorption due to the central wire. To this end we varied the length of wire 
inside the end of the tip and found the optimum value to be 180nm. The radius 
of the wire at the very end of the tip is 6.8nm. We then repeated the previous 
calculations, the results being shown in figures 3| to 



Comparing figures ^ to ^ to figures || to |7 
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shows the improvement in the 
magnitude of the transmission when using the new tip. For the 35nm aperture 
the improvement is about two orders of magnitude, falling to about one order of 
magnitude for the 25nm aperture. The general increase in signal strength has 
allowed us to perform the calculation for a smaller aperture than before, of only 
17nm radius. Also it is clear that the resolution for the coaxial tip is equal to if 
not slightly better than the conventional tip. The shape of the graphs at larger 
object radii is clearly much flatter than for the normal tip. We expect the fields 
to be strongly conflned in the region between tip and sample. The more tightly 
the fields are confined, the smaller the spot on the surface illuminated and so 
the less effect adding dielectric material to the object at large radii has on the 
transmission. This would seem to suggest fields more tightly confined for the 
coaxial tip compared to the normal one. 

Figure |ll] shows the transmission coefficient as a function of aperture radius 
for both types of tip and shows that for the range of apertures from about 65nm 
to 20nm the coaxial design offers significant improvement. 

5 Conclusions and future work 

In this paper we have presented the first attempt to apply the transfer matrix 
method to calculate the properties of the near-field optical microscope. We 
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have calculated the transmission coefficient for a range of aperture radii as well 
as estimating the resolution by comparing the transmission data for slightly 
different objects under the tip. 

Whilst our calculations may lack the sophistication of others that have been 
performed in this field, ours are the first, to the best of our knowledge, to use 
an adaptive mesh scheme to follow accurately the shape of the tip. There is no 
reason, in principle, why this co-ordinate transformation cannot be applied to 
other methods. Many other techniques such as Girard's method or the Finite 
difference time domain method rely on a real space discretisation at some point 
and so could directly benefit from the adaptive mesh we have developed here. In 
the future we also hope to remove the constraint of working within cylindrical 
symmetry and so calculate images of objects directly. 

We have also presented the first calculations for a near-field microscope with 
a coaxial tip. These results clearly show the potential for such a tip to enhance 
the magnitude of the transmitted signal even at very small aperture sizes. The 
only theoretical limit on the size of such a tip is the need for the central wire to 
have a diameter larger than the skin depth. 

Experimentally, some progress has been made towards realising such a tip. 



Fischer etal. \ 19 investigated the properties of a conical glass tip manufactured 
with a silver wire taper placed at the centre. However, difficulties in obtaining 
an unbroken metal core prevented him from using the tip for imaging. More 
recently, he has suggested a tetrahedral tip coated with metal except for 
one edge as an alternative which may offer the advantages of the coaxial design 
but without the manufacturing difficulties. It is hoped that calculations of the 
kind outlined here can be of use in analysing the merits of alternative SNOM 
designs such as these. 
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Figure 1: The tip-surface scattering system is contained within a cylinder. 
Waves can be incident from top or bottom, but the boundary conditions at 
the side must be specified. 
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Figure 2: A typical mesh used to represent the SNOM tip. 
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Figure 3: The mesh shaded to indicate the presence of the metal coating the 
tip (Diamonds), the substrate (circles) and the object being studied (squares) 
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Figure 4: Disc or ring objects (ring shown) can be placed under the tip. 
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Figure 5: Comparison of transmission coefficients for ring and disc objects. 
Aperture radius 48nm 
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Figure 6: Comparison of transmission coefficients for ring and disc objects. 
Aperture radius 35nm 
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Figure 7: Comparison of transmission coefficients for ring and disc objects. 
Aperture radius 25nm 
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Figure 8: Comparison of transmission coefBcients for ring and disc objects. New 
Coaxial tip of aperture 35nm 
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Figure 9: Comparison of transmission coefBcients for ring and disc objects. New 
Coaxial tip aperture 25nni 
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Figure 10: Comparison of transmission coefficients for ring and disc objects. 
New Coaxial tip aperture 17nm 
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Figure 11: Transmission coefficient against aperture radius for the conventional 
and coaxial tips. The inset is a plot of the log of the transmission against radius 
in the critical region 
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